Geometrical effects on energy transfer in disordered open quantum systems 
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We explore various design principles for efficient excitation energy transport in complex quantum sys- 
tems. We investigate energy transfer efficiency in randomly disordered geometries consisting of up to 20 
chromophores to explore spatial and spectral properties of small natural/artificial Light-Harvesting Complexes 
(LHC). We find significant statistical correlations among highly efficient random structures with respect to 
ground state properties, excitonic energy gaps, multichromophoric spatial connectivity, and path strengths. 
These correlations can even exist beyond the optimal regime of environment- assisted quantum transport. For 
random configurations embedded in spatial dimensions of 30 A and 50 A, we observe that the transport effi- 
ciency saturates to its maximum value if the systems contain 7 and 14 chromophores respectively. Remark- 
ably, these optimum values coincide with the number of chlorophylls in (Fenna-Matthews-Olson) FMO protein 
complex and LHC II monomers, respectively, suggesting a potential natural optimization with respect to chro- 
mophoric density. 



PACS numbers: 

Inspired by recent observations of Environment- Assisted 
Quantum Transport (ENAQT) in biological light-harvesting 
systems iT-TOl, one problem of fundamental and practical rel- 
evance is to engineer excitonic energy migration in disordered 
materials and nano-structures by exploiting the interplay of 
quantum effects and environmental interactions. Specifically, 
one might develop radically novel design principles by ma- 
nipulating delocalized exciton dynamics through engineering 
coherent couplings of free Hamiltonian and/or the environ- 
mental interactions to generate optimal efficiency. If one can 
fully understand the fundamental microscopic processes in- 
volved IlL JL2], it will be in principle possible to use structure 
for steering the ultrafast migration of excitons via quantum in- 
terference effects. Thus the excitons could be guided around 
defects and through interfaces in disordered systems with po- 
tentially important application in photodetection, bio-sensing, 
and photovoltaic light-harvesting LI 3 J . 

In this work, we investigate underlying geometrical and 
dynamical physical principles for efficient energy transport 
in generic random multichromophoric structures typically far 
from symmetrized configurations considered in Refs llflllTSl . 
Our work is motivated based on the recent studies on disor- 
der biological light-harvesting systems such as Fenna- 
Matthews-Olson (FMO) of green sulfur bacteria |1, 8 |, reac- 
tion center (RC) of purple bacteria 1 16], and light-harvesting 
complex II of higher plants 1 17 1. In particular, we investigate 
generic disordered configurations of up to 20 chromophores. 
We examine the energy transfer efficiency of the uniform 
distributions of the random light-harvesting complexes, as 
well as important subsets of these samples, encapsulated in 
a sphere of fixed diameter. For spatial sizes of 30A and 50A, 
corresponding to FMO and LHCII dimensions, we find trans- 
port efficiency increases as a function of the number of chro- 
mophores and saturates at optimal values of 7 and 14 sites 
respectively. Remarkably, these are indeed the number of 



chlorophylls for FMO and LHCII monomers implying a po- 
tential optimization with respect to chromophoric density. 

Moreover, we compute the energy transfer efficiency (ETE) 
of the uniform distributions of the random complexes consist- 
ing of 7 chromophores encapsulated in spherical dimension 
of various diameters interacting with a phononic bath. This 
allows us to compare our results with the performance of a 
well-characterized natural LHC such as the FMO complex. A 
similar study has been performed in Ref. fTSi, which have 
concludes that the FMO geometry with such high efficiency 
is extremely rare. In contrast to conclusions drawn in the 
Ref. [18], here we find that FMO performance is not rare 
if one uses a more accurate dynamical model, an appropri- 
ate measure of efficiency, and limits the search space to those 
configurations with similar compactness level as the FMO it- 
self. We use ETE as yield function for performance of these 
random systems and use a time-convolution master equation 
(TC2) 1 19 1 to estimated energy transport efficiency beyond 
perturbative and Markovian regimes ||9l|T0]|. We extensively 
explore the effect of spatial compactness on the transport ef- 
ficiency and robustness due to its potential significance as re- 
cently reported in ifTOl . Our simulations show that for larger 
size complexes the optimal configurations are not robust with 
respect to angular orientation of dipole moments in contrast to 
highly dense systems such as FMO. 

To explore potential geometrical patterns, we examine pos- 
sible spatial and energetic correlations among 10^ random 
multichromophoric samples embedded in a fixed diameter 
sphere ranging from 30 A to 100 A. We find that significant 
statistical correlations can exist among very high- and low- 
efficient samples even in suboptimal regime of chromophoric 
density, beyond the robust ENAQT regime |TOl. Specifically, 
we find certain structural similarities among various Frenkel 
exciton Hamiltonians of the high/low efficient random config- 
urations with respect to ground state properties, excitonic en- 
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ergy gap structures, and spatial connections. Moreover, we in- 
troduce a new measure of spatial path strength that can quan- 
tify an underlying structural mechanism for the performance 
of light-harvesting systems in a given dimension. 

This paper is organized as follows: In section II, we de- 
scribe our physical model and review the non-local master 
equation that was derived and analyzed in Ref. (9] to effi- 
ciently calculate ETE in multichromophoric systems. In sec- 
tion III, we examine the variation of ETE with the number 
of chromophores in fixed dimensions. In the section IV, we 
explore the role of chromophoric density in energy transfer 
as a dominating parameter [ilOj . In the subsequent sections 
we explore the roles of other physical parameters in transport 
including, ground state energy properties, average excitonic 
energy, structure and strength of chromophoric connectivity. 



phonon bath. The overall dynamics of the system is given by 
Liouvillian equation: 



dpsjt) 
dt 



-ih{[Htotah PSB{t)])ph = {^total[PSB{t)])ph 



(2) 

where psB denotes the system-bath state, and {...)ph indicates 
an average over phonon bath degrees of freedom. The Liou- 
villian superoperator Ctotai is the sum of superoperators £5, 
Cph and Cs-ph associated to Hs, Hph and Hs-ph respec- 
tively. The time evolved density operator of the multichro- 
mophoric systems in the interaction picture can be expressed 
by the propagator: 



pW = ('7+exp / Ctotai{s)ds )php{^) 



(3) 



I. THEORETICAL MODEL OF RANDOM 
MULTICHROMOPHORIC SYSTEMS 

Based on the generalized Bloch-Redfield master equations 
introduced by Cao |20|, we have recently derived the well- 
known time nonlocal master equation TC2 without making 
the usual weak system-bath coupling assumption |9|. By pro- 
viding an error analysis, we could show that TC2 can be em- 
ployed for highly efficient while reliable estimation of en- 
ergy transfer efficiency in light-harvesting complexes for both 
weak and intermediate system-bath coupling strengths and 
memory time scales. Here, we summarize the main steps of 
our approach. More technical details can be found in Ref. [9 1 . 

The dynamics of a photo synthetic system interacting with 
surrounding scaffold protein and solvent can be understood 
by starting from a general time evolution formulation of open 
quantum systems. The total system-bath Hamiltonian can be 
expressed as 



total 



= Hs + Hph + Hs- 



ph 



(1) 



where 



Hs 
Hph 
Hs-ph 



The phonon bath is modeled as a set of harmonic oscillators. 
Here \j) denotes an excitation state in a chromophore spatially 
located at site j. The diagonal site energies are denoted by e^s 
that include reorganization energy shifts \j = hj0^dP- j2 
due to interactions with a phonon bath; dj^^ is the dimen- 
sionless displacement of the (j, ^)th phonon mode from its 
equilibrium configuration. The strengths of dipole-dipole in- 
teractions between chromophores in different sites are rep- 
resented by Jjk- The operators Si = and Bj = 
— huj^dj^^Qj^^ are system and bath operators. Here, we 
assume that each site is linearly interacting with a separate 



where O denotes the interaction picture representation for an 
operator O. If we expand the above time-ordered exponential 
function, we arrive at the Dyson expansion for time evolu- 
tion of the density operator. This expansion involves multi- 
time bath correlation functions {B{ti^)...B{ti^^))ph, for any 
B = Bj. According to generalized Wick's theorem, for a 
system interacting with a bosonic (photonic and/or phononic) 
bath these higher order bath correlation functions can be ex- 
actly described by various combinations of only two-time cor- 
relation functions of the form Yl^ ^{Xj^B(ti^)B(ti^)) , where 
X+ is the index ordering operator. Assuming such Gaussian 
property for bath fluctuations, the most general approach to 
solve the master equation ^ is to utilize a path integral for- 
malism, leading to HEOM |21 1. However, such general ap- 
proach is impractical for our purpose as the computational re- 
sources required for simulating the energy transfer dynamics 
of photosynthetic complexes grow significantly with increas- 
ing size of the system, and with decreasing bath cutoff fre- 
quency and ambient temperature. 

We derived a numerically efficient method for simulation 
of complex excitonic systems by incorporating some physi- 
cal approximations in addition to Gaussian property to map 
quantum dynamics into a single solvable time-nonlocal equa- 
tion, see Ref. |9|. The main approximation is involved with a 
special truncation of higher-order bath correlation functions. 
Specifically, we assumed that generalized Wick's expansion 

{B{U,)...B(ti.^))ph = l[i^+mjB{U,)) (4) 



pairs l,k 



can be approximated as 



{B{U,)...B{U,J) ^ {X^B{t,)B{t2)){I^B{tk,)..B{tk,J)(5) 

for ti > t2 > ... > t2n- In other words, we disregard some 
fast decaying terms in the generalized Wick's expansion, but 
keep the slow decaying leading terms such that a two-point 
correlation can be factored out. This approximation can be 
understood phenomenologically by noting that two-point cor- 
relation functions Cj{t — ti) = {Bj (t)Bj (t 1 ) ) typically decay 
over a characteristic time, e.g., for a Lorentzian spectral den- 
sity and at high temperature T, the correlation functions decay 
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FIG. 1: The average energy transfer efficiency of random configu- 
rations consisting of different numbers of chromophores. The green 
[blue] dots represent the average ETE for 10^ random multichro- 
mophoric complexes embedded in a sphere of diameter 30 A [50 A]. 
The ETE increases with the number of sites nearly saturating at 7 and 
14 which coincide with the number of chromophores in the FMO 
and LHCII monomers with spatial dimensions of 28. 5 A, and 47 A 
respectively. These results imply a potential optimization of chro- 
mophoric density in nature. 



exponentially as A(2//3 — ^7)6"^*^^"*^^ where 7"^ is the re- 
laxation time of phonons and (3 = {kBT)~^ < 7"^. The 
approximation ^ is valid in the limit of time scales t longer 
than 7~^. As we see below, in calculating ETE we are con- 
sidering an integration over time therefore this approximation 
can provide good results even if 7 is not large or in another 
word we are in a relatively non-Markovian regime. 

Using this approximation we arrive at the time non-local 
master equation, TC2, as: 

d 

-p{t) = Csp{t) + Ce-hP{t) (6) 

where £e-/i = - ^Ls{b')01' •}-^^rap{|^rap)(trap|, .}, 
h.c. stands for Hermitian conjugate, and {,} is the anti- 
commutator ( [, ] is the communtator). The vector \trap) 
denotes the state of the site (bacteriochlorophyll (BChl)) 
connected to the reaction center. The term Ce-h = 
-Ej^Ls{b')0'l'-} - rtrap{\trap){trap\,.} captures two 
different competing electron-hole pair recombination pro- 
cesses that determine the energy transfer efficiency of light 
harvesting complexes. The first process, loss, is due to dissi- 
pation to the environment at each site that happens within the 
time- scale of 1 ns. This adverse environmental effect guaran- 
tees that the energy transfer efficiency has a value less than 
one. The second recombination process, trap, is due to suc- 
cessful trapping at one or more reaction center(s). 

A biologically relevant function for exploring the perfor- 
mance of light-harvesting complexes is the ETE as defined in 
Refs. 10^3^1, that is the total probability of exciton being 
successfully trapped: 



r] = 2rtrap {trap\p(t)\trap)dt (7) 
Jo 

We provide a formal derivation of the energy transfer effi- 
ciency in Ref. [9 |. ETE measures the likelihood of successful 
trapping, weighted by trapping rate: it quantifies the excita- 
tion availability whenever the reaction center is ready to op- 
erate within a period much shorter than the exciton life-time 
II22II . Note that this definition is very different than the first 
passage time |[l8l|24l, which quantifies the time-scale of first 
arrival of the exciton to the trapping sites. The latter definition 
is not necessarily correlated with the efficiency of quantum 
transport. In other words, for certain quantum processes the 
first passage time can be very short, compare to all other time- 
scales, but the transport mechanism could be still inefficient. 
In such cases, the excitations are typically delocalized over the 
regions that have very small overlaps with the reaction center 
and thus dissipate into the environment. 

The primary motivation for the truncation of correlation 
functions introduced above is to lead us to a special time non- 
local master equation that is solvable in the frequency domain. 
To estimate the regimes of the applicability of this method, 
one would ideally need to account quantitatively for the errors 
introduced by the generalized Wick's expansion truncation. 
We presented an approximate estimation of such inaccuracy 
in Ref. |9| for computing energy transfer efficiency by defin- 
ing an upper bound for the error as 

At] = 2rtrap\ / {trap\p{t) - pTC2{t)\trap)dt\ (8) 
Jo 

where p{t) is the exact density matrix of the system and 
pTC2 (t) is the solution to the TC2 master equation. Note that 
an exact account of errors in various regimes of interest is 
equivalent to calculation of the general evolution of the den- 
sity operator of the system that we intend to avoid. In Ref. O, 
using a combination of phenomenological and analytical ap- 
proaches, we found an approximate error bound for weak and 
intermediate system-bath couplings and bath memory time- 
scales, thus quantifying the reliability of our approach in such 
regimes. We also tested the accuracy of our method for by 
examining its predication in simulating quantum dynamics of 
FMO complex at room temperature compare to HEOM as a 
general benchmark ||25]| . We showed that oscillatory time evo- 
lution of the population of BChls in the FMO using our ap- 
proach are relatively close to those predicted by HEOM [9 |. 

Here, we would like to apply TC2 master equation to ex- 
plore the interplay of structurally-induced quantum coherence 
and environmental interactions for a large number of small- 
size random light-harvesting complexes with different chro- 
mophoric density, we choose the magnitudes of dipole mo- 
ments similar to the FMO, but with arbitrary random orienta- 
tions, site energy, and locations bounded in any given diame- 
ter of the spherical space. Our results hold for other materi- 
als with different dipole moment magnitudes by an appropri- 
ate renormalization of distances in spherical coordinate. We 
set the nearest neighbor distances by a lower bound of bA 
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FIG. 2: The effect of chromophoric compactness on ETE for different system-bath coupHng strength. We simulate ETE for 10^ different 
configurations of 7 chromophores embedded in a sphere of fixed diameter d ranging from 30 A to 100 A. The initial and target chromophores 
are located on the north and south poles of the spheres. The ETE histogram is depicted for reorganization energies of A = Ocm~^ (blue), 
A = 35cm~^ (green) and A = 350cm~^ (red). At A = 35cm~^ and d = 60 a uniform distribution in all classes of ETE is observed. 
Below d = 60 the samples are mostly high efficient and above d = 60 they become mostly low efficient. For A = and 350, this transition 
happens around d = 40 indicating that optimality and robustness of ETE can both be enhanced by the appropriate level of environmental 
fluctuations. The FMO complex belong to an ultra-high performing population represented by green bars in d = 30 implying the highest 
degree of efficiency and fault-tolerance in such chromophoric density and noise level 



due to intrinsic limitation of the dipole-dipole approximation 
|26 |. As we would like to explore the generic behavior of 
multichromophoric systems, here we do not account for spa- 
tial constraints due to particular size of each chromophore or 
a specific scaffold protein. We assume that the initial excita- 
tion and trapping sites are located at the surface of the sphere 
encapsulating a given configuration. That enables us to avoid 
a large amount of trivial optimal solutions in the configura- 
tion spaces; i.e. those in which the original donor and final 
acceptor sites are nearest neighbors. Thus the excitation has 
to travel through the entire length of the multichromophoric 
complex and generally experience multi-path quantum inter- 
ference in the regime of interest with intermediate system- 
bath coupling strength. 



A. Optimal number of chromophores 

Here, we explore the dependency of ETE on the number of 
chromophores for small light-harvesting complexes. In Fig. 
[T] we plot ETE for random complexes consist of 2 to 20 chro- 
mophores embedded in spheres of diameter d = 30A, and 
d = 50 A with environmental parameters given in fTl\. The 
average ETE for 1000 random configurations is computed for 
fixed diameters and number of chromophores. The standard 
deviation in these samples is depicted in the Fig. [T] inset. We 
observe that ETE increases monotonically by increasing the 
number of chromophores for different compactness levels. It 
is remarkable that for d = 30, which is the same as the di- 
ameter of the FMO complex, the number seven represents the 
minimal set of chromophores necessary to obtain high effi- 
ciency of 98%. Although, slightly higher ETE can be obtained 
by additional sites, that would be practically inefficient con- 
sidering the amount of work required to form such extremely 
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FIG. 3: Three samples of chromophoric arrangements from the 10^ random configuration in Fig. [2]for d = 60A and A = 35cm~^. These 
samples have distinct energy transfer efficiencies including very low efficiency of less than 1% in (a), intermediate efficiency of about 50% in 
(b), and high efficiency of around 93% in (c). It can be observed that simple geometrical considerations of spatial coordinates of chromophoric 
dipole moments cannot fully account for the significant discrepancies in their efficiencies. 



dense multichromphoric system for a marginal improvement 
in ETE. Here, the minimalist nature of natural selection might 
be at work: complexity is added until high efficiency and ro- 
bustness is attained, and no further. A similar behavior can 
be seen for d = 50 which coincides with the spatial size of 
LHCII. In such a distance, the ETE also saturates by increas- 
ing the number of sites and reaches to its optimal value for 14 
chromophores. Any extra chromophore may improve the av- 
erage efficiency by less than 1%. Ironically 14 is the number 
of chromophores of the LHCII monomers in higher plants. 

We expect that such saturation of ETE with respect to the 
number of mutlichromophores to occur for other range of 
compactness levels. It is of great importance if similar com- 
parisons with larger natural photosynthetic complexes can be 
demonstrated. That would imply a potential natural optimiza- 
tion with respect to the number of chromophores participat- 
ing in quantum transport. Similar studies for larger artifi- 
cial light-harvesting complexes could be of significant value 
for estimating the minimal number of chromophores needed 
to achieve a desired efficiency considering the physical and 
chemical spatial constraints in a realistic environmental con- 
dition. Such studies are beyond the scope of the current 
manuscript and will be undertaken in subsequent works. Next, 
we explore the role of chromophoric density by considering a 
fixed number of molecules in various spatial dimension. 



II. FUNDAMENTAL ROLE OF CHROMOPHORIC 
DENSITY 

We investigate the efficiency of random light-harvesting 
complexes, sampled from uniform distributions, embedded in 
a sphere of given diameter, d, ranging from 30 A to lOOA. For 
each compactness level defined by a fixed diameter, we cate- 
gorize the population of 10^ random arrangements in various 
classes based on their respective ETE and reorganization en- 
ergies. Figure [2] shows the histograms of such populations for 
various multichromophoric diameters for three different val- 
ues of reorganization energy chosen from the three different 



regions of system-bath couplings strength including, interme- 
diate regime, e.g., A = 3bcm~^ (green bars), fully coherent 
regime with virtually no environment 0cm ~^ (blue bars), and 
strong environmental interactions, e.g. 350cm~^ (red bars). 

Let's first examine the results associated to those complexes 
living in the similar environment as FMO (green bars in Fig. 
|2]). Observe that an overwhelming amount of random con- 
figurations have efficiencies comparable to the FMO complex 
for compactness level of about d = 30 A. Thus, in contrast 
to conclusions drawn in the Ref. |18|, the FMO performance 
is not rare if one limit the parameter space to those configu- 
rations that have a similar compactness as the FMO protein 
complex. Also, here we are using ETE as a measure of per- 
formance which capture long-time behavior on the same order 
as trapping time- scale. In addition, we use a TC2 dynamical 
equations that can go beyond Haken-Strobel model and cap- 
tures relaxations as well as pure dephasing process. 

Our results in Fig. [2] introduce an additional degree of the 
robustness of the FMO. It implies excellent tolerance with re- 
spect to the BChls locations, provided that its boundaries and 
environment are not changing radically. This robustness is dif- 
ferent than relative insensitivity of FMO transport efficiency 
to dipole orientations and site energies that were reported by 
us in Ref. |10|. Notably, as we increase the size of ran- 
dom complexes from 30 A to 100 A, for fixed A = 3bcm~^ 
(green bars), the histograms are drastically changing from be- 
ing sharply picked at high ETE to have sharp spectrum at low 
ETE, with average ETE dropping monotonically from 94% 
in d = 30A to 3% in d = lOOA. In the intermediate sizes, 
d = 50 A to d = 70 A, we observe that these samples are more 
or less evenly distributed in all efficiency levels with mean 
value of about 50 % efficiency for d = 60 A. This implies a 
smooth transition in ETE standard deviation, as a function of 
an effective parameter fi'^/d^, with optimal values within the 
range ofd = 60 A. 

One might expect that simple geometrical patterns can fully 
describe the variations of ETE in a given spatial dimension. 
For example, it is intuitively expected that straight line ar- 
rangements of the chromophores from initial excitation to 
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FIG. 4: The average (a) and the standard deviation (b) of ETE in the samples given in Fig. |2] The average ETE drops monotonically 
by increasing the diameter of the sphere from d = 30A to d = lOOA. The standard deviation shows maximum around d = 60A (40) 
for A = 35cm~^ (0 or 350). In such compactness regimes, the diverse populations of random configurations from very low to very high 
efficiency implies lack of robustness due to significant involvement of other physical parameters in guiding the exciton migration ETE, beyond 
the dominating factor of chromophoric density. 




Diameter 

FIG. 5: The overlap between the Hamiltonian ground state and the 
trapping site 3 is presented for the ratio of the average top 100 ef- 
ficient random samples over bottom 100 configurations. Although 
high efficiency geometries have larger overlap, especially at dimen- 
sion d = 60A, the overlap of ground at trapping site cannot fully 
account for the large discrepancies among top 100 efficient and 100 
low efficient configurations. 



trapping site should lead to more efficient configurations in 
a fixed dimension. However, by inspection of the actual loca- 
tions of the chromophores for three random samples associ- 
ated to diameter of d = 60 A we note that these apparently 
similar samples have significantly different energy transfer 
efficiencies, see Fig. [3] Thus, other structural and dynami- 
cal correlations should play important roles in discriminating 
among various samples with respect to ETE. We consider a 
variety of possible scenarios in the following sections. One 
potential important factor is the effect of environmental inter- 
actions on the shape of the ETE histograms and its transition 
from mostly efficient to mostly inefficient regimes. 

To study the impact of ENAQT on the above phenomenon, 
we simulate ETE for random chromophoric configurations in 
two extreme environmental cases in Fig. [2] Blue bars rep- 
resent the ETE for the ideal case of isolated systems and red 
bars show ETE for the systems that are strongly interacting 
with their surroundings (e.g, one order of magnitude stronger 
reorganization energies). The general features of the ETE his- 



togram persist but transitions occur in smaller dimensions for 
both of these cases around d = 3bA to d = 45A. Thus 
the existence of ETE statistical transition is independent of 
ENAQT. This phenomenon is essentially a direct manifes- 
tation of quantum dynamics driven by the internal Hamilto- 
nian, but its effect is modulated by reorganization energy. The 
ENAQT phenomenon can be seen here by noting that those 
configurations operating at (A = 35cm~^), represented by 
green bars, always have higher ETE at all compactness levels. 

Overall, by careful inspection of these results two main 
questions arise: How does the ETE behave as a function of the 
chromophoric density? What are the possible classical and/or 
quantum correlations, in the spatial and energetic structure of 
these random multichromophoric geometries, discriminating 
ultra high or low efficiencies in any fixed diameter? We ad- 
dressed the former question in details in the Ref. 1 10] by ex- 
amining the variation of the average transport efficiency for 
random configurations of up to 20 chromophores in two dif- 
ferent compactness level. We observed that the chromophoric 
density of the FMO complex and LHCII to be around the op- 
timal values for spatial dimensions ofd = 30 A and d = bOA 
respectively. Here, we investigate the latter question. 

In the following sections, we explore underlying structural 
and physical principle(s) for very high- or low-efficient 7- 
chromophoric configurations in any of the histograms in Fig. 
|2] beyond the dominating factor of compactness. In the first 
step of our analysis of Fig. [2j we compute the average ETE 
over all random samples in various spatial dimensions to em- 
phasize that very weak or strong environments are suboptimal 
in all ranges, see Fig. |4](a). In Fig. |4](b) the standard devi- 
ation of the average ETE is plotted for the same compactness 
levels highlighting the diversity of configurations with a max- 
imum pick at the range d = 40 io d = 60A for different 
reorganization energy values. 

III. GROUND STATE ENERGY OVERLAPS 

Based on the definition of ETE one generally expects that 
the overlap of lowest exciton state (ground state) with trapping 
site to be a good indicator of any potential correlations among 
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FIG. 6: The histogram of excitonic energy levels for 100 efficient samples and 100 lowest samples chosen from all 10^ configurations 
embedded in spheres with different fixed diameters ranging from d — 30 A io d — 100 A in Fig. [2] It can be observed that for the entire range 
of spatial dimensions low efficiency samples have energy gaps significantly larger than the phonons average energy 63 cm~^, and thus cannot 
use the bath as a energy sink to enhance exciton funneling in a fixed excitation limit. On the contrary, for top 100 highly efficient samples, we 
see that their average energy gap is sharply picked around 100 cm~^ or less for dimensions larger than d — 40 A. We note that for d — 30 A 
there is significant insensitivity with respect to these energy mismatches. 



highly efficient samples. This overlap is indeed very high for 
the FMO complex, about 0.94. In order to test this hypothesis 
we have diagonalized the free Hamiltonian of the top 100 high 
efficiency and bottom 100 low efficiency random samples, and 
calculated the average overlap between ground excitonic state 
with trapping. Fig. [5] shows the ratio of such average overlaps 
for highest 100 samples to lowest 100 samples. We note that 
some global correlations exist as the ground state overlap is 
bigger by a factor up to 4.5 for the high efficient samples in 
most compactness levels with maximum enhancement around 
d = 6OA. However, this correlation is incomplete, since it 
cannot account for huge discrepancies in ETE within each 
distribution. In particular, this measure cannot accurately dis- 
criminate ultra high (low) efficient samples (top/low 10 con- 
figurations) from the rest of geometries. Moreover, it cannot 
directly account for the role of environment. 



IV. EXCITONIC-PHONONIC ENERGY CONVERGENCE 

One measure that can potentially capture the role of envi- 
ronment in discriminating among various configurations is the 
compatibility of exciton energy gaps with phonon energies. In 
the weak system-bath couplings the multiphonon transitions 
are not common and most of environment induced dynamics 
is driven by single phonon transitions. In such regimes, those 
configurations that have exciton energy gaps comparable with 
the energies of phonons can easily use the bath as an energy 
sink. This can lead to enhancement of funneling toward trap- 
ping sites by passing extra energy to phonons in a fixed exci- 
tation manifold. Using bosonic distribution function, a rough 
estimation of the average energy of a phonon can be simply 
evaluated by ignoring the protein chemical energy, as follows: 



/Q^c^cjJ(cj)cj/(exp(/3cj) - 1) 
i^dujJ{uj)/{eMM-^) 
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FIG. 7: The excitonic energy levels of high- and low-efficient sam- 
ples. The green (blue) line shows the average energy levels of top 
(bottom) 100 samples. The green (blue) dashed-line shows the en- 
ergy levels standard deviation of top (bottom) 100 samples. The en- 
ergy of phonons at T = 298 °K and 7 = 50 cm~^ has an average 
of 64 cm~^. It can be seen that the high-efficiency samples have 
energy levels in the range of phonon energies, therefore facilitating 
energy exchange with the bath and enhancing the ENAQT process. 



where J{uj) is the Lorentzian spectral density. Using the 
above relation the average energy of a single phonon is about 
64 cm~^, assuming bath with cutoff frequency of 50 cm~^ 
at room temperature. We also compute the average energy 
gaps, g, for the bottom 100 and top 100 efficiency samples. 
The distribution of such samples based on their average ex- 
citon gaps in various compactness levels is shown in Fig. [6] 
We observe that bottom 100 samples (grey bars) have average 
energy gaps between 200 to 400 cm~^ with large standard de- 
viation of about 200 cm~^, see Fig. [t] Thus, there is typically 
an energy mismatch between typical exciton gaps and phonon 
energies. On the contrary, the top 100 efficient samples have 
energy gaps of about 100 cm~^ with smaller standard devia- 
tion of about 50cm~^. Consequently, there is a considerable 
chance of system-bath energy exchange facilitating exciton 
transfer. A more accurate description of the exciton band gap 
correlations should include multiphonon transitions for strong 
system-bath couplings that can be captured in Forster theory 
12811 or modified Redfield theory 1 29] . However, our data sug- 
gest that one- or two- phonon transitions that are relevant in 
the intermediate regimes can be captured by our dynamical 
equations Eq. |6] 



V. SPATIAL CONNECTIVITY 

In the classical regime of incoherent hopping, a simple de- 
sign principle for the multichromophoric systems, coupled via 
dipole-dipole interactions, to achieve high ETE would be to 
align chromophores on a straight line connecting the initial 
and target sites. In Fig. [3] we observe that relatively similar 
patterns might lead to configurations with very different ETE. 
By enlarging the spherical size encapsulating the samples, the 
average inter-chromophoric couplings become weaker, there- 
fore the quantum coherence would be attenuated. Thus, we 
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FIG. 8: The average distance of the 5 intermediating chromophores 
from the Z axis connecting the initial and target chromophores for to- 
tal number of random configurations as well as top/bottom high/low 
efficient samples. The average distances for various ETE classes be- 
come more distinct by increasing the radius of the spheres. For larger 
spheres, the attenuation of quantum effects enforces the high ETE 
configurations to have narrower structures of intermediating chro- 
mophores. The average distance of 9.8 A for 100 high-efficient sam- 
ples at radius 50 A versus 17 A for the total samples indicates this 
fact. This behavior is not significant for more compact geometries or 
low efficient samples 



expect that the measure of proximity to the axis connecting 
the initial and target chromophores, z, takes lower (higher) 
values for high (low) efficiency samples. To this end, we ex- 
amine this measure defined as the average distance of the 5 
intermediating chromophores from the Z axis. This numerical 
study, plotted in Fig. [8] confirms such intuition to some ex- 
tent. The average distance of total number of random configu- 
rations from the Z axis is plotted, as well as such distances for 
top 100 (high-efficient) and bottom 100 (low-efficient) sam- 
ples. It can be seen that the average distances for these differ- 
ent classes become more distinct by increasing the dimension 
of the sphere. For diameter 100 A, the top 100 samples on av- 
erage are closer to Z axis compare to the total populations with 
a ratio of 0.6. On the other hand, this feature is less important 
for more compact structures and for low-efficient configura- 
tions suggesting that simple geometrical consideration per se 
cannot fully determine the performance of an excitonic energy 
transfer system. 



VI. PATH STRENGTHS 

One important structural feature of light harvesting com- 
plexes is the spatial connectivity of their constituent chro- 
mophores. Here we demonstrate that such geometrical pa- 
rameters in the site base could provide an underlying physical 
explanation for the vast diversity of high or low energy trans- 
fer efficiency of random configurations embedded in a fixed 
volume. Fig. [2] We first define the concept of spatial path and 
its path strength between initial and trapping sites. 

Generally, for any m-chromophoric system, there are 
5^fe^o^ spatial paths each including different cou- 
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FIG. 9: The histogram of the strength of 326 possible paths connecting initial and trapping states for the very top-/low- efficiency configuration 
(chosen from 10^ random 7-chromophoric light-harvesting complexes in a fixed diameter 30 A to 100 A presented in Fig. [2]). In each fixed d 
the right (left) histogram demonstrates the number of paths for a given range of path strength, in energy units of cm~^, for the most (least) 
efficient sample. A considerable energy gap of about 1500cm~^ in the left panel is observed that is one order of magnitude bigger than those 
in the right panel for all volumes. There are a few strong paths that dominate the energy transport in the low efficient samples. In contrast, a 
large number of paths are contributing to quantum transport for high performing samples in the right panel. These alternative paths could help 
avoiding quantum localization due to static or dynamical disorders leading to a substantial enhancement of energy transport efficiency. 



pling combinations of0<k<m — 2 chromophores connect- 
ing initial and target sites. Label the initial site as 1 and the tar- 
get site as 7. A path defined by interconnecting chromophores 
{ci, c/c|l < Cj < 7}, has path strength, hc^^,,,^ck^ defined 
as the inverse of the total time scale of going through the path 
from the chromophore 1 to 7. This time scale h~^^_ ^_^^^ is given 
by the sum of the inverse of coupling strengths between neigh- 
boring sites given by the off-diagonal elements of the Hamil- 
tonian H 

k-l 

KL.c, = \Hi,cA-' + E l^c,,c,+J-' + \H,,,7\-' (10) 

and the trivial path of sites 1 and 7 direct connection has 
strength Hi j. For a light harvesting complex consisting of 
7 chromophores, there are 326 spatial paths each including a 
different coupling combination of one to five chromophores. 
In order to study any potential relationship between the path 
strength (as defined above) and ETE, we compute the path 
strength for all the paths for the most and least efficient con- 
figurations, from 10000 random 7-chromophoric samples in 
a given sphere with diameter 30 A to 100 A. The results are 
demonstrated in Fig. [9j where in each fixed d the right (left) 
histogram demonstrates the number of paths for a given range 
of path strength in energy unites of cm~^ for the most (least) 
efficient configuration. 

We note that there are significant path strength gaps of 
about 1500cm~^ in the left panel (ultra low ETE samples) 
which are one order of magnitude bigger than those in the 



right panel (ultra high ETE samples with gaps of less than 
150cm~^) at all compactness levels. Incidently, for low effi- 
ciency random configurations there are a few paths that com- 
pletely dominate the energy transport whereas the high effi- 
cient samples contain a large number of, more or less, uni- 
formly distributing spatial paths that contribute to quantum 
transport. One reasonable explanation is that light harvest- 
ing complexes with a few effective paths for exciton trans- 
port are more susceptible to large energy mismatches due to 
static or dynamical disorders leading to quantum localization. 
By contrast, those complexes with a large number of active 
spatial path become more robust to disorders or defects by 
providing alternative routes of transport. Moreover, the multi- 
path systems can better exploit the environmental fluctuations 
to overcome energy mismatches which otherwise lead to ex- 
citon localizations. A similar robustness to defects due to 
redundant paths has been recently observed experimentally 
for energy transfer in certain artificial light-harvesting com- 
plexes ||30l[3Tl and numerically simulated using an incoherent 
Forster model. These artificial systems are synthesized from 
self assembly of tobacco mosaic virus scaffold protein in disc 
and rod geometries which in the latter up to thousands of chro- 
mophores can be positioned spirally. Overall, the concept of 
path strength provides a rather straight forward explanation 
for the diversity of configurations in the ETE histograms pre- 
sented in the Fig. [2] 

To investigate this phenomenon beyond the highest or low- 
est efficient configurations, we have performed statistical 
studies on different ensembles of the samples, see Fig. [9] 
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FIG. 10: (a) The number of pathways with path strength larger than lOOOcm"^ statistically averaged over m highest ETE samples. The 
horizontal axis shows the number of high ETE samples included in the statistical computation. It can be observed that top 500 high efficient 
configurations have no dominant pathways. If we enlarge the number of configurations beyond top 500, i.e., including samples with interme- 
diate ETE, a few paths with strong strength will appear, (b) The number of spatial pathways with strength larger than 1 000cm ~^ statistically 
averaged over m lowest ETE samples. The horizontal axis shows the number of low ETE samples contributed to the statistical average. 
In sharp contrast to histogram (a), the energy transfer in the low efficient configurations is dominated by a few channels with very strong 
path strength. This property is attenuated when we increase the number of statistical samples to include more efficient ETE configurations. 
These results demonstrate that the path strength, defined by Eq. [T0| is a suitable measure for quantifying the geometrical mechanism for the 
performance of light-harvesting systems in a given spatial dimension. 



(a) and (b). We observe that for the low efficient configura- 
tions the excitonic energy is carried through a few dominant 
pathways whose strength is one or two order of magnitude 
larger than the rest of pathways. Thus, for m top or bot- 
tom samples we count the number of the paths with strength 
larger than lOOOcm"^ as a measure for having a few dom- 



inant paths. In Figs. 10 (a) and (b) we compute the aver- 
age of this measure for m high and low ETE samples (m = 
10 - 90, 100 - 900, 1000 - 5000) and for d = 30, 60, 90A. 
Remarkably, it can be observed that all of top 500 high effi- 
cient samples do not have any number of dominant paths, see 
Fig. [To] (a). However, as we compute the average over a 
larger number of high efficient samples some stronger paths 
will emerge. As expected, the inverse behavior is observed 



for low efficient samples presented in Fig. 10 (b), where the 
low ETE configurations have some dominant pathways that 
disappear as we enlarge number of low ETE samples. These 
results demonstrate that the path strength, defined by Eq. \T0 
capture the spatial correlations among top/low ETE structures 
and thus can provide an underlying geometrical description 
for the efficiency of light-harvesting complexes. 



VII. CONCLUSION 

We studied distributions of random arrangements of chro- 
mophores in volumes of various diameters as a function of 
chromophoric density, reorganization energy, and their in- 
terplay. We demonstrated that the chromophoric density 
play a major role in determining the performance of few- 
chromophoric systems embedded in various spherical vol- 
umes. We found that for random mutlichromophoric systems 
with fixed sizes 30 A and 50 A, 7 and 14 chromophores are 
indeed the minimum numbers to achieve a high ETE, similar 



to those values for FMO and LHCII respectively. After our 
original observations of these optimal chromophoric numbers 
in biological LHC ifTOl , a similar result has been recently re- 
ported 1321 . 



We observed significant statistical correlations in the low- 
and high-end efficient random structures with respect to aver- 
age exciton energy gaps and multichromophoric spatial con- 
nectivity at the intermediate system-bath couplings. More- 
over, we have investigated possible statistical correlations be- 
tween quantum entanglement and performance of random 
chromophoric configurations. Specifically, we simulated the 
dynamics of excitons for many high/low ETE samples and 
calculated the entanglement measure based on Ref. |33|. Our 
simulations showed no particular patterns from which one can 
conclude that the presence of quantum entanglement is a sig- 
nificant indicator of efficient energy transport. 



Overall, we observed two distinct parameter regimes for 
efficient energy transport: Ideally one can design molecular 
configurations to reside in the so-called Quantum Goldilocks 
regime (TOl |34l where there is an energy- scale convergence 
for coherent and incoherent processes leading to optimality 
and robustness. However, if such regime is practically unac- 
cessible, here we found alternative scenarios for potential op- 
timal material design beyond Goldilocks regime where certain 
configurations could still have high energy efficiencies, but 
they become very sensitive to defects or environmental fluc- 
tuations. These model studies could be useful for design of 
structured molecular aggregates such as self-assembled dyes 
on tubular J-aggregates 1 35 , 36] and virus-based templates ag- 
gregates with potential applications to photo- 
voltaic devices, photosensing, and biological sensing. 
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